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Abstract 

We fix the long-standing ambiguity in the 1-loop contribution to the mass of a 1+1-dimen- 
sional supersymmetric soliton by adopting a set of boundary conditions which follow from 
the symmetries of the action and which depend only on the topology of the sector considered, 
and by invoking a physical principle that ought to hold generally in quantum field theories 
with a topological sector: for vanishing mass and other dimensionful constants, the vacuum 
energies in the trivial and topological sectors have to become equal. In the two-dimensional 
N — 1 supersymmetric case we find a result which for the supersymmetric sine-Gordon model 
agrees with the known exact solution of the S-matrix but seems to violate the BPS bound. 
We analyze the nontrivial relation between the quantum soliton mass and the quantum BPS 
bound and find a resolution. For N = 2 supersymmetric theories, there are no one-loop 
corrections to the soliton mass and to the central charge (and also no ambiguities) so that 
the BPS bound is always saturated. Beyond 1-loop there are no ambiguities in any theory, 
which we explicitly check by a 2-loop calculation in the sine-Gordon model. 
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1 Introduction 



The calculation of quantum corrections to the mass of solitonic objects has been a sub- 
ject of intense interest in the past |IJ and has recently been revived in the light of a re- 
cent breakthrough in our understanding of non-perturbative dynamics in 3+1-dimensional 
super symmetric (susy) gauge theories 0. One of the most important ingredients of the 
non-perturbative analysis of such theories is the duality between extended solitonic objects, 
such as monopoles or dyons, and point-like particles. Another important ingredient is the 
concept of the BPS spectrum — the particles whose masses are proportional to their charges. 
Due to the supersymmetry algebra BPS states may be annihilated by the action of some 
of the supersymmetry generators and hence give rise to a smaller number of superpartners 
("multiplet shortening"). Therefore the BPS value of the mass becomes a qualitative, rather 
than just a quantitative property M. 

Susy models in 1+1 dimensions provide a valuable nontrivial testing ground for these 
concepts. Solitons in these models are examples of BPS states whose mass at the classical 
level is proportional to the topological charge. The calculation of quantum corrections to 
the masses of these states has been subject to a long controversy. A number of one-loop 
calculations have been performed yielding different contradicting results ||, ||, |6|, [7|, |8|, |9|, [10], 



11, 12, 13 



The one-loop corrections to the mass of a soliton are given by 

m$ = I E - - F ) - \ E - - F ) + SM (1) 

where u B,F are the energies of the small bosonic (fermionic) fluctuations about the classical 
soliton solution, Q are the corresponding energies of the linearized theory in the trivial vac- 
uum, while 5M is the counterterm which can be obtained from the expression for the classical 
mass of the soliton in terms of unrenormalized parameters by expanding into renormalized 
ones | 15| . In order to make the above sums well defined, spatial boundaries are temporarily 
introduced to make the entire spectrum discrete. 

One can identify two sources of ambiguities. As was discussed in ||, imposing different 
spatial boundary conditions on the small quantum fluctuations around the classical soliton 
gives different, sometimes even ultraviolet divergent, results. In this paper we present an 
analysis which answers the question: which boundary conditions are to be used in the one- 
loop calculation? The answer to this question can be found if one re-examines the original 
formulation of the problem. We consider the vacuum energy as a functional of the boundary 
conditions. We then single out a class of boundary conditions which do not introduce surface 
effects — the topological boundary conditions. They close the system on itself. There is 
a trivial (periodic) as well as a topologically non-trivial (with a Moebius-like twist) way 
of doing this. This definition of boundary conditions does not rely on a semiclassical loop 
expansion. We do not separate the classical part of the field from its quantum fluctuations; 
rather, the boundary conditions are imposed on the whole field. One then infers the correct 
boundary conditions for the classical part as well as for the quantum fluctuations from this 
single general condition. 

Another source of ambiguity, as was pointed out recently fl4|], is the choice of the ul- 
traviolet regularization scheme. The dependence on the choice of regularization scheme can 
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be understood as a peculiar property of those quantities, such as the soliton mass, which 
involve a comparison between two sectors with different boundary conditions, i.e., different 
topology. Indeed, the difference of the vacuum energies in the two sectors measures the mass 
of the soliton. The one-loop correction is then given by a sum over zero-point frequencies 
in the soliton sector which is quadratically divergent. A similar sum in another (the trivial) 
sector is to be subtracted in order to get an expression which is finite if written in terms 
of the renormalized parameters. It turns out that due to the bad ultraviolet behavior of 
both sums the result depends on the choice of the cut-off |L4 ] . Cutting both sums off at the 
same energy in both sectors J7], ||, |IU| or taking equal numbers of modes in both sectors 
fl6| , ||, |l|, [13| leads to different results. To add to the confusion, some authors do not include 
bound states and/or zero modes when they consider equal numbers of states in both sectors.^ 
In this paper we propose a simple way of reducing the ultraviolet divergence of the sums over 
the zero-point energies, which eliminates the sensitivity to the ultraviolet cutoff. Instead of 
calculating the sums we calculate their derivative with respect to the physical mass scale 
in the theory. The constant of integration can be fixed by using the following observation: 
the vacuum energy should not depend on the topology when the mass is zero. This is the 
physical principle which allows us to perform the calculation unambiguously. It should be 
viewed as a renormalization condition. 

From a practical point of view we need to use our condition that the vacuum energy 
functional does not depend on topology at zero mass only in the one-loop calculation. How- 
ever, to preserve the spirit of our approach, we formulate this condition, as we do with our 
boundary conditions, for the full theory regardless of the semiclassical loop expansion. The 
mass that needs to be taken to zero is the physical, renormalized, mass scale. From this 
point of view this condition is a trivial consequence of dimensional analysis if we work with a 
renormalizable theory where all the physical masses are proportional to one mass scale. All 
the masses, including the soliton mass, vanish then at the same point, which is the conformal, 
or critical, point in the theory. Another way to look at our condition to fix the integration 
constant is to consider the Euclidean version of the 1+1 theory as a classical statistical field 
theory. Then the mass of the soliton is the interface tension between two phases. As is well 
known [17 the interface tension vanishes at the critical point; moreover, it vanishes with the 
same exponent as the inverse correlation length. 

In section |2| we demonstrate our new unambiguous method of calculation using as an 
example the bosonic kink. We show that, as argued in [14|, the correct result corresponds to 
mode number cutoff. (The same conclusion was recently reached for nontopological solitons 
in 3+1 dimensions |I8| .) In section ^| we apply our analysis of the topological boundary 
conditions to the case of an N = 1 susy soliton, where it leads to nontrivial consequences. 
We analyze the relation of our results to the BPS bound in section |j. In section |5| we analyze 
the N = 2 susy solitons and conclude that the one-loop corrections vanish completely. In 
section ^] we redo the 2-loop corrections for the case of the bosonic sine-Gordon soliton 
|19| , p0[1 , paying this time close attention to possible ambiguities, and find that no ultraviolet 
ambiguities appear. The ultraviolet ambiguity is thus purely a one-loop effect which leads 
to the interesting conjecture that it may be formulated in terms of a topological quantum 

1 See for example the textbook Q, eq. (5.60). Actually in this reference the result of Jl6| is obtained but 
this requires neglecting a boundary term in the partial integration of (5.63). 
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anomaly. 



2 Eliminating the one-loop ultraviolet ambiguity using 
a physical principle 

In this section we present a general analysis regarding the calculation of the soliton mass 
which will help us eliminate the ultraviolet ambiguity discussed in [f[4l| . We consider the 4 
theory (kink) as an example, but the arguments can be applied to the sine-Gordon theory 
as well. The crucial property of these models from which our boundary conditions follow is 
the Z 2 symmetry — > — 0. 

Let us take a step back from the actual calculation and try to define the mass of the 
soliton before we do the semiclassical expansion. We start from the observation that the 
soliton carries a conserved charge — the topological charge. This means that we can define 
the mass of the soliton as the difference between the energy of the system with nontrivial 
topology and the energy of the system with trivial topology. This definition coincides with 
the definition based on path integrals in the topological sector which are normalized by path 
integrals in the trivial sector ||19|| . The topological charge of the system is determined by 



the conditions at the spatial boundary. We view the vacuum energy as a functional of the 
boundary conditions. In general, a boundary condition could induce surface effects associated 
with the interaction of the system with the external forces responsible for the given boundary 
condition. We would like to avoid these contributions. There is a class of boundary conditions 
which do not produce such effects. These are what we call topological boundary conditions, 
which identify the degrees of freedom at different points on the boundary modulo a symmetry 
transformation. In our case there are two such possibilities: periodic and antiperiodic. These 
are dictated by the internal Z(2) symmetry: <ft(x) — > (—l) p <ft(x), p — 0,1. Crossing the 
boundary is associated with a change of variables leaving the action invariant: </>(— L/2) = 
(— l) p 0(L/2). The system behaves continuously across the boundary, only our description 
changes. In effect such boundary conditions do not introduce a boundary, rather they close 
a system in a way similar to the Moebius strip .0 

We spent so much time on this, perhaps, trivial point in order to make the choice of 
boundary conditions for the theory with fermions clear. The analysis of fermions, however, 
will be postponed until the next section. In the literature a large number of other boundary 
conditions have been considered, both in the bosonic and in the fermionic sectors but from 
our perspective they all introduce surface effects or are even inconsistent. 

It should now be clear that the mass of the soliton can be defined as the difference of the 
vacuum energy with antiperiodic and with periodic boundary conditions when the volume 
L — > oo. This definition does not rely on the semiclassical expansion. Returning now to 
the semiclassical calculation we see that at the classical level the equations of motion select 
the trivial or the soliton vacuum configuration depending on the topology. Less trivially, 
we see that at the one-loop level the boundary conditions that should be used for the small 
fluctuations about the soliton configuration should be antiperiodic. We must point out that 



2 Though topological boundary conditions do not induce boundary effects, finite volume effects vanishing 
as L — > co are of course present. Such effects will be discussed in section 0. 
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the choice of boundary conditions does not affect the result of the calculation in the purely 
bosonic case. We shall nevertheless use antiperiodic boundary conditions in the soliton sector 
in this section to be faithful to our nonperturbative definition of the soliton mass. We shall 
see in the next section that for fermions the choice of the boundary conditions becomes 
crucial. 

A few points about the classical antiperiodic soliton should be stressed here. The topo- 
logical boundary condition reads 

<f>(-L/2) = (-in(L/2) and <f/(-L/2) = (-1)^/2), (2) 

where the nontrivial sector is selected when p — 1. Note that the derivative with respect to 
x, 4>', must also be antiperiodic in the soliton sector. The classical soliton solution <p{x) can 
be viewed as a trajectory of a particle with coordinate <j) moving in time x in the potential 
—V(<f>). The particle is oscillating about the origin = with a period which depends on 
the amplitude. When the period is equal to 2L the trajectory during half of the period is the 
antiperiodic soliton satisfying the boundary conditions (0). The endpoints of the trajectory 
need not necessarily be the turning points. For example, the particle at time —L/2 can start 
downhill at some <f>o with nonzero velocity, then pass the point at the same height on the 
opposite side, i.e., — <fio, going uphill, then turn and after that at time L/2 pass the point — </>o 
again, but going downhill. Clearly, for this trajectory, @ is satisfied, whereas restricting the 
usual soliton solution centered at x = to the interval (—L/2, L/2) would lead to a solution 
for which <p(-L/2) = -<j>{L/2), but <f>'(-L/2) = +<f>'(L/2). When L -> oo the turning points 
of the trajectory come infinitesimally close to the minima of V{(f) and we recover the usual 
L = oo soliton. 

Next, we want to address the problem of the ultraviolet ambiguity in the one- loop cal- 
culation. To summarize the beginning of this section 

M = E X - E , (3) 

where E p , p = 0, 1 are the energies of the system with the periodic and antiperiodic boundary 
conditions of @ respectively. At the classical level this gives M c \ ~ m 3 /(3A) for the kink, 
where m is the mass of the elementary boson at tree level, and A the dimensionful coupling 
constant. The order K correction is due to the fact that boundary conditions change the 
spectrum of zero-point fluctuations. Two factors are responsible for the ultraviolet ambiguity 



discussed in [14]. One is the fact that the terms in the sums over zero-point energies grow 
making each sum strongly (quadratically) ultraviolet cutoff dependent. Second is that one 
has to compare the spectrum in two different vacua. Taking all the modes below a certain 
energy in both systems leads to a different result than taking equal numbers of modes. It 
would be nice if there was a parameter in the theory whose variation would continuously 
interpolate between the two vacua. This is not possible due to the topological nature of the 
difference between the vacua. However, we can identify a certain value of the dimensionful 
parameters of the theory at which the vacuum energy should become independent of the 
topology. This will be one ingredient of our calculation. 

Another ingredient is the observation that one can reduce the ultraviolet divergence of 
the sums of zero-point energies by differentiating w.r.t. m. The terms in the differentiated 
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sums become then decreasing and as a result the sums (now only logarithmically divergent) 
can be unambiguously calculated. But the price is that we need to supply the value for the 
integration constant to recover the function from its derivative. This can be done using a 
physical principle that relates the energies of the two vacua at some value of the mass. One 
must realize that the difference in the energies arises because of the nontrivial potential for 
the scalar field. If this potential vanishes the energies of the two vacua become equal. In the 
absence of the potential the mass m of the boson is zero and the soliton disappears. Therefore 
the constant of integration over m is fixed by the condition that the energy difference between 
the two vacua must vanish when m — > 0. A subtlety here is that m should be sufficiently 
large compared to 1 /L so that finite volume effects can be neglected. The limit m — > should 
be understood in the sense that the mass approaches 0(1/ L), where L is large. Then the 
difference between the vacuum energies must be 0(1/ L). Also note, that other dimensionful 
parameters in the theory should be scaled accordingly when m — > 0, e.g., A/m 2 = const in 
the \<j) 4 theory. 

We want to relate the mass of the soliton to other parameters of the theory. The relation 
to the bare parameters mo and A will contain infinities. The infinities in the relation of phys- 
ical quantities to the bare parameters in this theory should be eliminated if we renormalize 
the mass mo = m + 5m, where 

3A r°° dk 1 

Sm = / . (4) 

2m J-oc 2tt + m 2 V ; 

With this renormalization of mo tadpole diagrams vanish. In the 4 theory the physical pole 
mass of the meson differs from m by a finite amount — y / 3A/(4m) ||14|| ; however, it suffices 
to use m for our purposes. If we rerun this analysis for the sine-Gordon theory, the tadpole 
renormalized mass m would, at one-loop order, coincide with the physical meson mass. If 
we use this renormalized mass in the expression for the soliton mass we get an additional 
one-loop counterterm, 8M, 



where 



ml 


m 


3A 


= 3A H 


m 2 5m 


3m r 


A 


~Y J- 



If, (5) 

,,. ,„,. f°° dk 1 
SM= — - — = — / . 6 

° ' 2tt y/k 2 + m 2 V ; 

Now we differentiate the well-known expression for the one-loop correction to the 
soliton mass with respect to the mass m 

dM^ _ d(5M) + ly^duj n ly, dCu n 
dm dm 2 „ dm 2 ~^ dm 

For the spectrum u n in the trivial sector one obtains 

k„L = 2rm. (8) 



dm Jk 2 + m 2 ' 
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For the soliton sector 



dm 



ki + m 2 



m + k n 



dm 



m + 



1 kl 



S\k n ] 



kl+m 2 \ Lm 

k n L + 5{k n ) = 2im + n, 



(9) 



where we used the fact that S(k) depends on m only through k/m to convert the derivative 
w.r.t. m into the derivative w.r.t. k. 

We convert the sums over the spectrum in (0) into integrals over k using 



COO rjh 



and 



£/(*») 



(10) 



(ii) 



These expressions follow from the Euler-Maclaurin formula which is valid for a smooth 
function f{k) vanishing at k = oo. In our case /(fc) = dui/dm and /(&) = du/dm satisfy 
these conditions. From the Euler-Maclaurin formula one can also see that in the naive 
calculation with f(k) = u the ambiguous contribution, which comes from regions 8/L at the 
ultraviolet ends of the integration interval, is non-vanishing due to the fact that f(k) grows 
with k. 



\ \ \ w 
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Figure 1: The left- and the right-hand sides of the equation 5(h) = 27m + tc — kL are plotted 
schematically by solid lines in the case of the 4 kink (two bound states). The dashed 
line represents the value of S(k) without the discontinuity 2Tie{k). Observe that with this 
discontinuity the mode numbers n — — 1, should be left out, while the spectrum of allowed 
values of k is not affected. 



We can use the following expression for the phase shifts 5(k) in the case of the kink in 
b A theory: 

5(k)= ^27r-2arctan J^'^ W), ( i2 ) 
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where we added the term 2irs(k) to ensure that 8(\k\ — > oo) — > 0. As a result, 8(k) is 
discontinuous at k = 0: 8(G±) = ±2n. It is then easy to see (Fig. [I]) that for n = — 1 
and n = the equation ([)]) does not have solutions for k. It is pleasing to observe that 
this defect is matched by the existence of two discrete modes: ujq = (the translational 
zero mode) and u-i = y/3m/2 (a genuine bound state). To these discrete modes we can 
assign (somewhat arbitrarily) those "unclaimed" n's. That this matching is not a coinci- 
dence follows from Levinson's theorem: 8(0±) = ±7rnd s , where n^s is the number of discrete 
solutions. Since the discontinuity in 8(k) is an integer multiple of 2tt, it does not change 
the spectrum of the allowed values of k (Fig. [I]). This spectrum near the origin is given by 
kL = . . . , — 3ir, —it, 7T, 3tt, ... up to 0(1/ L) for 5 either with or without the 2ire(k) term in 
(f[2f), and the values of 8'(k) and f(k) on this set of fc's are also not affected. 
Putting now all the pieces together we obtain 

dm dm 2^ dm 2m J-oo2ix K J K J 

This formula is universal. Substituting for the 4 theory the particular values for 8M from 
(|) and 8(k) from QTJ) we findfl 

dM^ 3 m 2 fQO dk 1 ^3 

+ 



dm 2 J-oo 2tt (A; 2 + m 2 ) 3 /2 4 

3m 2 dA; 1 13 



2 7-oo 2tt VA; 2 + m 2 (m 2 + 4P) Ay/S 2vr' ^ 

Integrating over m and using that = when m = we obtain the result for the one-loop 
correction to the kink mass which was previously obtained using mode number cutoff EBL [14| 



M«=m . (15) 

3 Fermions, supersymmetry, and topological boundary 
conditions 

In this section we shall extend the ideas introduced in the previous section to theories with 
fermions, and in particular theories with supersymmetry. The following analysis can be 
applied to any N = (1, 1) supersymmetric theory with Lagrangian 

£ = ~\d^<P - l -U 2 (<P) - ~ (m + U'(<j>)W) , (16) 

where U(4>) is a symmetric function, admitting a classical soliton solution. For example, for 
the kink U{<j>) = J\/2(<p 2 - ^). We use {7^, Y} = ^ u with g 00 = -1, g 11 = 1, and ip is 
3 Use /^(l + x 2 )- 1 /2(i + 4x 2 )~ 1 dx = 2>/37r/9 and J^(l + x 2 )- 3 / 2 ^ = 2. 
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a Majorana spinor: ip = ipU^ = ip T C with = —{ r y^) T C. The action is invariant under 
5(f> = ei[) and 8ijj = (fi<f> — U)e. 

First of all we want to identify the class of topological boundary conditions. The discrete 
transformation taking <ft — ► —0 must be accompanied by ip — > 73^ with 73 = 7°7 1 to leave 
the action invariant. From this symmetry transformation we obtain topological boundary 
conditions 



<K-L/2) = (-in(L/2), 

H-L/2) 



//(-L/2) = (-1)V(V2), 
( 73 )^(L/2), p = 0,l 



gives a trivial periodic vacuum while p 



(17) 

1 selects a nontrivial soliton 



The value p 
vacuum. 

As one can see, the reasons behind our choice of the topological boundary conditions do 
not include supersymmetry. The same arguments apply to any theory with a Yukawa-like 
interaction between fermions and bosons. From this point of view it is very gratifying to 
discover that the p = 1 topological boundary condition fllTD preserves half of the supersym- 
metry of the Lagrangian (|T6|). An easy way to see that is to consider the Noether current 
corresponding to the supersymmetry 



Integrating the conservation equation d^J^ = over space we find 



at ~ atJ- L /2 ax 1 



J\x) 



L/2 
-L/2 



(19) 



where the r.h.s. is simply the total current flowing into the system. Using the boundary 
condition (1171) we obtain 



J\x) 



L/2 
-L/2 



+ f/) 7 V - {-H + f/)7 1 7 3 ^ 



L/2 



;i+73)(^>+^hV 



L/2 



(20) 



We see that the (1 — 73) projection of the supercharge Q is conserved. Note that different 
projections of Q are classically conserved on the soliton or the antisoliton background. The 
soliton with </>' + U = and ip — is invariant under a susy transformation with a parameter 
e if (1 + 7 x )e = 0. This means the projection P+Q (with P± — 1 ± 71) of the supercharge 
vanishes on the soliton configuration to linear order in the quantum fields. For the antisoliton 
P-Q has this property. This should be expected since the topological boundary condition 
does not distinguish between the soliton and the antisoliton. 

Similarly one can see that the topological boundary condition does not break translational 
invariance. The conservation equation for the stress tensor reads d^T^ = 0, where 



4 

In general, the non-conservation of total momentum is again due to the boundary term 

dP d f L / 2 



(21) 



at at J -l/2 



T 



01 



T 



11 



L/2 
-L/2 



(22) 
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We see that the defining property of the topological boundary condition, that it relates 
the fields at —L/2 and L/2 by a transformation leaving £ invariant, ensures momentum 
conservation. 

Note also that there is another Z(2) symmetry in the Lagrangian (|16|): ip — > (—l) q ip. 
This can be used to extend the set of topological boundary conditions (|T7|) to 

4>{-L/2) = (-in(L/2), 4>'{-L/2) = (-l)V(^/2), 

${-L/2) = (-m l3 ) p ^(L/2), p,q = 0, 1. (23) 

The values (p, q) = (0, 0) give a topologically trivial sector. The sector (0, 1) is also trivial, 
but the fermions have a twist. Two classically nontrivial vacua are obtained with p = 1 and 
q = 0,1. For p = 1 the two values of q correspond to the arbitrariness of the sign choice of the 
73 matrix, and are related to each other by space parity transformation ip(x, t) — > 7°-?/>(— x, t). 
Therefore one should expect E(l, 0) = £7(1, 1), which one can check is true at one-loop. 

As in the previous section we define the mass of the soliton as the difference of the energies 
E p of these vacua: M = E% — E Q . At the classical level one finds M = M c \, where M c \ is the 
classical soliton mass. The one-loop correction is determined by integrating 

dMW _ d(5M) l^du* 1 dug (1 ^ du* 1 du£\ 
dm dm 2 ^ dm 2 ^ dm \2 ^ dm 2^ dm) 1 ' 

over m. The expressions for the derivatives du^/dm and du^/dm are the same as in the 
bosonic case, see (H) and ([5]). In order to find the corresponding expressions for the fermionic 
frequencies we need to obtain the quantization condition for k n . For the trivial sector we 
have simply k n L = 2im. 

The nontrivial sector requires more careful analysis. The frequencies u are obtained by 
finding solutions of the equation 

(y^ + U') V> = (25) 

of the form ip = ip(x) exp{— iujt}. Multiplying this equation by {—^^d^ + U') we find 

(-d 2 + (U') 2 + 1 1 U"u)iP = 0, (26) 

where we used 4>' sol = —U((j) so \), which follows from the classical equation of motion for the 
soliton in the L — > oo limit. Projecting this equation using P± = (1 ± 7 1 )/2 we see that ip + 
(where ip± = P±ip) obeys the same equation as the bosonic small fluctuations, hence 

1p + = a+e -^+i(kx±8/2) when x _^ ±oc _ ^ 

The ip- component can then be obtained by acting with P + on (|25|) 

d ol °^ + (d 1 + U')^ + =0. (28) 
which together with (p7f ) gives the asymptotics of ip- 

^ = - 7 °q + fe J iT " e -faW(te±V2) = _^ a+e ±i6/2 e - M {k X ±5/2) when x _^ ±QO) (29) 

v /c ~\~ in 
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where 6(k) = — 2 arctan(m//c) and we used the fact that in the L — > oo limit U' — > ±m. 
Therefore the solutions of ([25]) have asymptotics 



V> = V>+ + V>- = (l - 7°e ±ie/2 ) a+e -*-*+i(**±*/2) when a; -> ±00, (30) 

where a + is the eigenvector of 7 1 with eigenvalue +1. 

Although one could continue the derivation without specifying the representation for the 7 
matrices (an exercise for the reader) we find it more concise to adopt a certain representation. 
The most convenient is the following one in terms of the Pauli matrices: 7 = — ir 2 , 7 1 = r 3 , 
and hence 73 = T\. It has two advantages. First, a + has now only an upper component. 
Second, the Majorana condition becomes simply if>* = ip and the equation (|25|) is real. In 
this representation (|3CiD becomes 




^ = * ae -i^+i{kx±5/2) when x ^ ±00j (3^ 



where a is a complex number. 

Now we impose the boundary condition ip(—L/2) = ^if)(L/2). The field if) in equation 
(E5|) must be real. This means that only the real part of ([J]]) need to satisfy the boundary 
condition. This condition should, however, be maintained for all t. Therefore, due to the 
oscillating phase exp(— hot), a complex equation must be satisfied 

e- l(kL+S) ( =r( len) (32) 



We introduced T = 73 in order to discuss briefly the following point. One could consider 
a more general boundary condition: if>(— L/2) = Tif)(L/2) with some matrix T (real in the 
Majorana representation we have chosen). One can see from (|32"D that certain boundary 
conditions cannot be satisfied, for example, the frequently employed [0, [H| |ll|, [12], |T3| 



periodic boundary conditions with F — 1. Equation ( p2|) provides an additional consistency 
check for our choice of boundary condition. 

With the topological boundary condition T = 73 = r 1; we find that equation ([32] ) is 
satisfied provided 

Q 

kL + 5 + - = 2vm + tt. (33) 
Using this quantization rule we find 

duj^ 1 / , dk 



m + kn-£ =-=_ m+ T e + (34) 




dm + m 2 V dm I ^k 2 + m 2 \ L m 



Now we convert the sums over modes into integrals. For the bosonic and fermionic sums 
in the trivial sector formula flit] ) applies and the sums cancel each other (no cosmological 
constant in the trivial susy vacuum). For the bosonic sum in the nontrivial sector we use 
again (pT[). For the fermionic sum a formula analogous to (|TTJ) applies with 5 + 8/2 instead 



of S, which follows from fl33|) . Again, as for the bosonic modes, due to the discontinuity in 5 



there are values of n which do not lead to a solution of (0). The remaining n lead to k 
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Figure 2: The left- and the right-hand sides of the equation 5 + 6/2 = 2im + tc — kL are 
plotted schematically in the case of the supersymmetric 4 kink (two bound states). The 
dashed line represent the value of 5(h) + 8(k)/2 without the discontinuity 27re(k). As in the 
bosonic spectrum the discontinuity leads to n = —1,0 mode numbers being skipped, while 
the spectrum of allowed values of k is not affected. 



values (namely, in the case of the kink: kL = . . . , —5n/2, — vr/2, 3tt/2, 7tt/2, . . ., see Fig. 
which in the limit L —>■ oo give a continuous integration measure Ldk/(2ir) near k = 0. 
Putting now everything into formula (pil) we find 



dM^ d(SM) 



00 dk 



dm dm 4m J-oo 2n 

The one-loop mass counterterm is given by [|14| 

m r°° dk 1 
27-00 2^^/W+ 



Vk 2 + m 2 6\k). 



6M 



m^ 



(35) 



(36) 



It follows from the renormalization counterterm 5m which is chosen to cancel the sum of 
the bosonic and fermionic tadpole diagrams. Substituting into ( pop we find dM^/dm = 
— l/(2n), hence 

M« = - — . (37) 
2tt 



This result differs from the one two of us have obtained previously ]I3J using a mode- 
number regularization scheme with the conventionally employed (but, as we have argued, 
untenable) periodic boundary conditions. 

In the special case of the supersymmetric sine-Gordon model, we can compare this result 
with the one obtained from the Yang-Baxter equation assuming the factorization of the 



S-matrix [21]. The mass spectrum is then given by |2! 

m n = 2Msin(n7/16), 



where 7 in the notation of ref. [22] is related to the bare coupling {3 through 

1 l-/3 2 /4vr 1 1 



7 



4/3 2 



167T 



(38) 



(39) 
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Expanding for n = 1 we find 

m 1 = ^ + 0( 7 3 ). (10) 



Since this is the lightest mass in the spectrum we identify it with the meson mass m = m\. 
Taking the ratio Mjm\ and using ( p9|) we obtain 

- = - + 0( 7 ) = -|-f + 0(/3 2 ). (41) 

The first term is the classical result, the second is the 1-loop correction. This means that 



the 1-loop correction to M following from the exact S-matrix factorization calculation p2 
is the same as our 1-loop result fl37|). 

The next question we address is whether such a negative correction is in agreement with 
the well-known BPS bound. This question has been subject to controversy and deserves a 
separate section. 



4 Quantum BPS bound, soliton mass, and finite size 
effects 

As was first realized by Olive and Witten ||, the naive supersymmetry algebra in a topo- 
logically nontrivial sector is modified by central charges. The susy generators for the N = 1 
model read in the representation of section 3 

Q± = P± [-@<P + Uh°^}dx = / + ± U)i^]dx, (42) 

J-L/2 J-L/2 

where P± is again (1 ± 7i)/2. Using canonical commutation relations we arrive at the 
following algebra 

{Q±,Q±} = 2H±2Z; {Q + ,Q_} = 2P, (43) 

where 

r L / 2 11 1 i 

h = j_ L j-^ + -(^) 2 + -c/ 2 +^^+V'-+V'-V';)-^ / V' + V'-]^ (44) 

P = [ L/2 W + k^' + + ^'_)}dx (45) 

J — L /2 — 
pL/2 

Z = / ip'Udx (46) 

J-L/2 

The central charge Z is clearly a boundary term. 

Let us comment on some subtleties in the derivation of (|43|). The Dirac delta functions in 
the equal-time canonical commutation relations can be written as 5(x, y) = J2mVm{x)Vm{y), 
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where rj m is a complete set of functions satisfying the boundary conditions of the correspond- 
ing field. For such S(x,y) one has: 

" « t)i(x _ m = ^/ 2 > * «-^> ,47) 

-L/2 2 

1/2 0(x)<S(x + L/2) = *(-*/2> **W2) (4g) 

-L/2 2 

For the bosons in the topological sector we need the — signs. For the fermions ip+{x) +ip-(x) 
one needs the + signs, but for the fermions i[>+(x) ~'4 ) -( x ) one needs the — signs. That some 
subtlety in the delta functions is present, is immediately clear if one considers the double 
integral / / dxdy f(x)d x 5(x — y)g(y), and either directly partially integrates the derivative 
d/dx, or first replaces d x 6(x — y) by — d y 5{x — y) and then partially integrates w.r.t. y. One 
gets the same result provided 

[L/2 rL/2 

fg{x)\xeB= dy f(x)S(x-y)g(y)\ XIE B+ dx f(x)5(x - y)g(y)\ yeB (49) 

J-L/2 J-L/2 

where the notation h(x)\ x£ B implies h(L/2) — h(—L/2). Naively, there is a factor of 2 missing 
in this equation, but with the more careful definitions of 5(x — y) for periodic or antiperiodic 
functions, consistency is obtained. 

With these delta functions one finds that the boundary terms in the {Q±, Q±} anticom- 
mutators reproduce the consistency condition Q4"9]), hence cancel, whereas in the {Q + ,Q_} 
anticommutator one finds the boundary term 

- % J dxip + (x)5 a (x - y)4>+(y)\ y eB - i J dxi)_{x)5 a {x - y)if)-(y)\ yeB (50) 

where 5 a (x — y) is the bosonic (antisymmetric) delta function defined above. These terms 
cancel if one uses our topological boundary conditions. In the {Q, Q} relations one does 
not encounter subtleties involving delta functions for the fermions because there are no 
derivatives of fermions in Q±. 

Since the operators Q± are hermitian, one finds that the following relation exists between 
the expectation values of operators H and Z: 

(s\H\s) > \(s\Z\s)\, (51) 

for any quantum state s. 

As we have already pointed out in section 3 Q2DD, only one linear combination of Q + 
and Q_ is conserved in the soliton sector: Q + ± Q_ for q — 1,0 respectively. Taking for 
definiteness q = we derive from (|4"2l) , fl4lf) the following commutation relations: 

i[H,Q + + QJ\ = 2[(<p + <p')ty + + iJ>_.) + U(i/> + -i/;-)] x=L/2 (52) 

t[P 1 Q + + QJ\ = 2[(0 + /)(^ + + V-)-^(V+-V'-)]«=£/2, (53) 

while the other linear combination Q + — Q_ commutes with both H and P.^ We also find 
that the operator Z does not commute with the Hamiltonian in the soliton sector: 

i[H,Z] = 2U ( p\ x=L/2 (54) 



4 In these relations one must use the proper definitions of the delta functions for the fermions since 
derivatives of ip± appear in H. As one can see, (p2) agrees with (EQ). 
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Let us examine carefully the meaning of this last result. Strictly speaking, it implies that 
Z is not a central charge. We shall show that this fact reflects a certain property of Z in 
a finite volume L. Let us ask the following question: what is the expectation value of Z in 
the soliton vacuum state | sol) ? It is clear that for any finite volume L: (sol|Z|sol) = 0, since 
neither a positive nor a negative value is distinguished by our boundary conditions. This 
is a consequence of the fact that H and Z do not commute. The ground state |sol) is an 
eigenvector of H, but it need not be and is not an eigenvector of Z. 

To make the next step clear it is convenient to use the following observation: The value 
of Z measures the position of the soliton. Indeed, for the classical configuration if the center 
of the soliton is at x = exactly then Z is M (up to 0(1/L)). If we now move the center 
of the classical soliton the value of Z will decrease and reach zero when the center of the 
soliton is at L/2. If we continue shifting the solution in the same direction, and bearing in 
mind its antiperiodicity, Z will become increasingly negative and it will reach the value — M 
when there is an antisoliton at x — 0. If we deal with a configuration which is a distortion 
of the soliton, the center of the soliton is not well defined but Z is and can give one an idea 
of where the soliton is. On the quantum level this corresponds to the fact that P and Z do 
not commute. If we act with P (this is our spatial shift) on an eigenstate of Z we generate 
a different eigenstate of Z and the expectation value of Z changes. 

Now, if we think of Z as (a nonlinear function of) the coordinate and P as the momentum 
of the soliton, the next step becomes clear. The vacuum is an eigenstate of P with eigenvalue 
0. Therefore, it is a superposition of the eigenstates of Z. Positive and negative values of Z 
enter with the same weight into this superposition (the corresponding eigenvectors can be 
obtained from each other by acting with exp(iLP) — a shift by L). Therefore (sol|Z|sol)=0 
for any finite L. 

It is, however, too early to conclude that the equation (|5T| ) does not lead to any condition 
on (sol | if | sol) apart from semipositivity. The expectation value of Z can be compared to an 
order parameter in a system with spontaneous symmetry breaking. It is only nonzero if the 
thermodynamic limit L — > oo is taken properly. We shall now analyze how the limit L — > oo 
must be taken in the case of the soliton. 

One can view the soliton as an almost classical particle (as long as the coupling constant 
is small) which is subject to Brownian motion due to quantum fluctuations. This is the 
meaning of the fact that H and Z do not commute: Z, or the position of the soliton, 
depends on time. If we wait for a sufficiently long time it will cover all possible positions 
and the expectation value of Z will be zero. However, it is obvious that for large L most 
of the time the soliton will spend away from the boundaries. It means that if one starts 
from a state of the soliton away from the boundary so that Z = M and limits the time of 
observation the expectation value of Z will remain close to M. 

How long does it take for the soliton to cover all the volume LI Since it is a random walk 
the distance from the original position is 0(yt) and it will stay away from the boundary 
if we restrict the interval t of the observation to t <C 0(L 2 ). Such a restriction will mean 
that we introduce an error of at most 0(1/L 2 ) in the energy due to uncertainty principle. 
This is small when L is large. Alternatively, one can do Euclidean rotation and consider 
classical statistical theory of an interface in 2 dimensions. The random walk in this case is 



the well-known roughening of the interface [17]. It leads to smearing of the interface to the 
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width of w = constyt/M, where M is the one-dimensional interface tension, which is the 
soliton mass. The correction to the tension turns out to be 0(l/t), as we have already seen. 
(More rigorously, the partition function of the wall is not just exp(— Mt) in the t — > oo limit 



but has a preexponent LyM/t due to the fluctuations of the interface. The factor L arises 
from the integration over the volume of the collective coordinate with a familiar measure 
a/ Mt, and an additional factor \jt comes from the determinant of nonzero soft vibrational 
modes of the interface p3|). 



Therefore the thermodynamic limit L — > oo in our system which gives nonzero (sol|Z|sol) 
corresponds to the energy measurement whose duration t is small compared to ML 2 . The 
bound (^) must apply to the result of such a measurement. The expectation value of Z in 



such a thermodynamic limit is what was calculated in [p| |14" 

d 



V(0soi(*))<?7 2 (*)> 



/oo d 
dx— 
-oo dx 

m f°° dk 1 . , . , 

= M cl m - - / — 2 = M cl (m). (55) 

The logarithmically divergent integral is exactly cancelled by the counterterm 5M ([36|) after 
the renormalization of m. 

The question we need to answer now is what is the value of (sol|iJ|sol)? We recall that 
our calculation of the soliton mass was aimed at finding the dependence of M, or, (sol|iJ |sol) 
on m, the renormalized mass. We were not able to determine a constant term in the unrenor- 
malized (sol|-ff|sol), but we knew that it must be subtracted to satisfy the renormalization 
condition M \ m =o = 0. To evaluate the l.h.s. of ( |5"TD we have to know this constant. In order 
to find it we must evaluate directly the sums of bosonic and fermionic frequencies in the 
soliton sector. Although the sums are quadratically divergent supersymmetry improves the 
situation. It requires that bosonic and fermionic modes come in pairs. Therefore we need 
to apply the Euler-Maclaurin formula to a function f[n) = (u>^ — w^)/2 which has much 
better behavior at large n. Using the spectral relations for the bosons (|9|) and the fermions 
(|33|) we find: 



<sol|tf|sol> =M cl (m ) + i E (^-^) = M cl (m ) + l rg^^v/FT 

= MJmo) + — VP + m 2 9(k) - - / — Vk 2 + m 2 9\k) 
87r —A 4 J -a 2ir 



A-^oo 



A m 1 r°° dk m , . 

Md(mo) + T - - n I oT /,,, i 2 - ( 56 ) 



4 2tt 2 J~oo 2w VF + 



m 



The constant A is the ultraviolet cutoff related to the number of modes: LA = 2tc(N+ 1/2) + 
0(1/ L), and we recall that 9 = —2 arctan(m / k) . The last integral is logarithmically divergent 



and is exactly cancelled by the counter-term 5M fl36|) when we renormalize mo = m + 5m. 
The same divergence appears in (sol|Z|sol) and is also removed when m is renormalized. On 
the other hand, the linear divergent term A/4 does not appear in (sol|Z|sol). In terms of 
the renormalized mass m and M c \{m) the left- and the right-hand sides of (^l|) are: 

A Tfl 

(sol|F|sol) = M cl + - - — and (sol|Z|sol) = M cl . (57) 
4 2tx 
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We see that the bound is observed by the soliton vacuum state with an infinite overkill due 
to a linearly divergent constant A/4. This constant is only nonzero in the soliton sector. In 
this sector the bound is not saturated, which shows that "what can happen — does happen" : 
the argument of Olive and Witten for the saturation of the bound is based on the "multiplet 
shortening" and does not apply to N — 1 susy solitons. 

In other words, to resolve the long-standing problem of the Bogomolnyi bound in the 
N — 1 supersymmetric soliton/kink model we must realize that the bound is imposed on 
the unrenormalized expectation value of the Hamiltonian (sol | H | sol). Both sides of equation 
(|5"T|) contain ultraviolet divergences. The supersymmetry ensures that quadratic divergences 
do not appear in (H). However, N = 1 is not enough and a linear divergence remains in 
the topologically nontrivial sector (there is also a logarithmic divergence, but it is matched 
on both sides, (H) and (Z)). This linearly divergent term is positive in accordance with 
the bound (|5T|). Note that this divergence is different from a cosmological constant (which 
vanishes because of supersymmetry) in that it is proportional to L° rather than L 1 . 

Since (sol | H | sol) is divergent even after standard renormalization of the mass m we need 
to use an additional renormalization condition to find the physical soliton mass. This is the 
condition M| m=0 = which we introduced. Therefore 

M = (sol|#|sol) - (sol|#|sol)| m=0 = (sol|#|sol) - j. (58) 

Our new renormalization condition is based on the physical requirement that the physical 
vacuum energy should not depend on topology in the conformal point m — 0. Therefore, 
in principle, it also requires a subtraction of the expectation value of H in the topologically 
trivial vacuum (0|if|0), or rather (0\H |0) — (0\H |0) | m= o- In the non-sup ersymmetric case this 
is essential to cancel the background bulk contributions linear in L, which are m-dependent, 
and which are the same in both sectors. However, supersymmetry ensures that (0|if|0) = 0. 
Finally, we find a negative finite quantum correction to the physical soliton mass given by 
(|3T|). However, it is not the physical mass M, but rather it is (sol | H | sol), to which the 
bound (|5l"D applies. As we shall see in the next section, in the case of N = 2 supersymmetry 
all corrections, even finite, are cancelled and M = (sol|if|sol) saturates the bound as the 
"multiplet shortening" demands. 

5 The N=2 case 

Consider the action for the following N=(2,2) susy model in 1+1 dimensions 

£ = -d^d^cp - W<9^ - U*U - ^(U'^p T ij°tp + U*'ip* T ij°tp*) (59) 

where <p and ip are complex, U = U(ip), U' = (d/d(p)U, and 7 1 = r 3 and 7 = — ir 2 with r 3 
and T2 the usual Pauli matrices. The action is invariant under 5ip = tip, 5(f* = ipe, Sip = 
pipe — U*e* and 5ip = —eftip* — Ue T irf with complex two-component spinors e, with e = eHj 
and x/j = ipU'-y . 

5 With 7 1 = T3 instead of 7 1 = t\ the diagonalization of the fermionic actions is easier. 



16 



In terms of the components ip = (f\ + if2)/V^ and ip T = the action reads 

1 



C 



For the superkink 



;<9 M y?i<9 / Vi 



1 



V\(if 2 



Mo/2A) 



£4ink 

while for super sine-Gordon theory 

^sine-Gordon = ?™o COs(\f\/2(j) / m Q ) 

In terms of y?i and y?2 the potential is given by 

A 



kink 



Z^O 2 



U*U R i 



sine— Gordon 



mo 
A 



COS if i + COSh <^ 2 

m m 



(60) 
(61) 

(62) 

(63) 
(64) 



This already shows that for the kink the trivial solutions if^ = ±jj, /\/X and the kink 
solution ip i = fx — (/<o/ vA) tanh/io^/V^ (and the antikink solution ip^ = —fx) of the 
bosonic kink model remain solutions of this susy model while f2 — 0. Because the potential 
is of the form (f 2 — [i 2 ) /2\)(f* 2 — Hq/2\) instead of (iff* — /io/2A) 2 , there is no U(l) symmetry 
acting on (if\, if 2) which can rotate the kink away. Hence there is a genuine soliton. 



In sine-Gordon theory, the trivial vacua are at ip^ = an d f^ 
for the vacuum in the topological sector we choose the solution 



(n + 2)71"'/// 



/y/\, while 



mj vA(4 arctan(exp 



mi 



7T 



This solution is antisymmetric in x, in agreement with the Z 2 topological boundary condi- 
tions of section 2.Q 

The transformation rules in component form read 



5f* = —iip* + e^ + iip*_e + Sip- 



if e + — ife_ — U e + 
—f e_ + fe + — U e_ 



(65) 
= 



and since for the soliton (f' sol + y/2U(ip so \/y/2) = 0, the solution if = (l/\/2)f soh = ip_ 
is only preserved by half the susy transformations, namely those with 1m e + and Re e_. 

In the sector without soliton we set if\ = (fx + rj and find then the following linearized 
field equations 



□77 - m r] = , Dif 2 - m f2 = 



ip'_ + m ip*_ 
ip' - m ipl 



1 ^+ - ^+ + ™oVh 

J -ip- — + 








(66) 



6 This solution is obtained from the one in |Q by the substitution ip = tp' + mir/\f\. Actually, the 
sine-Gordon model has Z symmetry, and we could choose solutions which interpolate between two other 
minima of the potential by constant shifts. 
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where m = V (4^ / V2) and we used that U(4^/V2) = 0. Hence all fields satisfy the 
same second-order field equation in the trivial sector, with a common mass m (for the kink 
m 2 = 2f4). 

In the sector with a soliton, we set (fx = (p so \(x) + rj, and find then the linearized field 
equations 

□77 - (U' 2 + UU") r] = ; ip + -il)'_ + U' ip*_ = 



<Psol 



Vsol 



Uip 2 -(U' 2 -UU") ^ = ; ^_ VI = (67) 



'/'sol 



Vsol 



Decomposing ^ + = Re^+ + ilm^p+ and similarly for , 0_, the fermionic field equations split 
into one pair of equations which couple Reip + to Reip-, and another pair which couple 
Im4>+ to Im^-- Iteration and the relation (p' sol = —\/2U lead to 

Re 4)1 -Re/p+- (U' 2 + UU") Re 4j + = , idem for Im ^_ 

Vsol 

Re 4)"_ -Re^_- (U' 2 - UU") Retp. = 0, idem for Im ip + (68) 

Hence, the real triplet rj, Re^ + and Imvp- satisfies the same field equation as the real scalar 
r) and the real upper component 4 ) + of the Majorana fermion in the N=(l,l) model, whereas 
the real triplet ip2, Irrupt and Re^p- satisfies the same field equation as the real spinor ?/>_ 
in the N=(l,l) model @. 



In principle one can directly determine the discrete spectrum of y?2 by solving the 
Schrodinger equation with V = (U') 2 - UU" f24|, but susy already gives the answer. In 
the N=(l,l) model, the spinors u± in 4>±( x ,t) — u±(x) exp —iut satisfy the coupled equa- 
tions (d x + U')u + + ituu- = and (d x — U')u- + icuu + = 0, where U = \/2U((p = (fi/y/2). 
Any solution for u + and u__ yields also a solution for f2 (with (f2 ~ w_), and any solution 
for if2 leads also to a solution for u + and w_ (with U- ~ (f2 and u + ~ (d x — U')(p2)- Hence, 
there are as many bound states for ^2 as for rj, namely <^2,b ~ (d x + U')i]b- (The zero mode 
u + ~ ip' sol does not lead to a corresponding solution for U- and y?2 since (d x + U')ip' soX = 0). 
The discrete spectrum of the fermions in the N=(2,2) model is then as follows: in the sec- 
tor with Retp+ and Reip- there are one discrete state with zero energy (Re^+ ~ <^soi) 
and bound states with energy u B (Reip + ~ Vb)- In the sector with Reip- and Imip- the 
same normalizable solutions are found. Thus, as expected, the massive bosonic spectrum 
of small oscillations around the soliton background is equal to the corresponding fermionic 
spectrum, and consists of massive quartets. There are also one bosonic (for translations) 
and two fermionic zero modes (zero energy solutions of the linearized field equations). The 
latter are proportional to the nonvanishing susy variations Sip ~ ^' so \ e i which are due to the 
susy parameters Re e + and Im e_ .|j The zero modes do not form a susy multiplet (there are 
two fermionic and only one bosonic zero mode) but this poses no problem as quantization 
of collective coordinates tells us that the translational zero mode does not correspond to a 
physical particle. The fermionic zero modes are due to translations in superspace, namely 
when the susy generators Q± act on the superfield $(x, 6) = ip so \. 



7 One can directly determine these fermionic zero modes from ( |67| ) by looking for time-independent nor- 
malizable solutions. One finds then Reip + ~ Imip- ~ exp[— J* U'((j> so i(x')/>/2)dx']. These functions are 
indeed proportional to <p'i. 
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The topological boundary conditions for the action (|5j|) are 

0(-L/2) = (-l)^(L/2); t//{-L/2) = (-iy<f>'(L/2); 

H-L/2) = (-l)«( 73 )^(L/2). (69) 

For the continuous spectrum with p — 1 and g = or q = 1 boundary conditions, the 
quantization conditions are 

7] : fcL + = 2nir + tt 
y?2 : kL + 5(k) + 6(k) = 2mr + tt 

Reij + ,Imip + ,Reip-,Imip_ : fcL + + = 2n7r + gvr (70) 

The one-loop corrections to the soliton mass, differentiated w.r.t. m, are then given by 
<9 m M (1) = d m 5M {1) 

+— / VF + m 2 (*'(*) + {5\k) + ^(fc)} - 2{<5'(fc) + -e\k)}) ^ (71) 

2m J-oo V 2 / 27r 

where the massive bound states do not contribute because they come in susy multiplets. 
The continuum states do not contribute either, because all phase shifts clearly cancel. The 
counterterm 5M vanishes in the N = 2 model since the rj and 02 tadpole give ^ and ^ times 
a fermionic tadpole, respectively. Note that all these cancelations are in fact independent of 
any particular regularization scheme since all integrands cancel. 

To decide whether the Bogomolnyi bound is saturated we now turn to the central charges. 
The super Poincare charges are obtained by the Noether method and read 



Q- = - [{Uif> + -il)*_<p-4) m + <p')dx (72) 



With topological boundary conditions the combination Q + — Q_ is conserved, but not Q+ + 

Q-, 

Q + + Q_ = -2 [U(iP+ - V-) + (r+ + r-)(<P + <p')] L/2 (73) 

There are no ordering ambiguities in these operators, and using equal-time canonical 
commutation relations one finds the following algebra for A± = Q + ± {Q+)* and B± = 
Q- ± Q*_ 

{A ± ,A ± } = ±2H-(Z + Z*); {A + , A_} = -(Z - Z*) 

{B±,B ± } = ±2H + 2ReZ; {B + ,B_} = Z - Z* 

{A ± ,B ± } = ±2P;{A + ,B_} = {A_,B + } = (74) 

The generators which are produced on the right-hand side are 

H = f°° ((pip* + (p'(p*' + UU* + iip+ip+ + iiplip-) dx 

r+oo , , 

P = lip' p* + ipip*' + + iifj*_ip'_ ) dx 

J —oo ^ ' 
roc 

Z = 2 U^'dx (75) 
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Since an N = 1 massless multiplet in D = (3, 1) (which is always without central charge) 
becomes a massive N = 2 multiplet in D = (1,1) whose (mass) 2 is equal to the square 
of the central charges, while the N=(l,l) susy algebra in D = (3, 1) becomes an N=(2,2) 
susy algebra in D = (1,1) with two central charges (the generators P2 and P3), it is clear 
why massive multiplets of D — (1,1) N = (2,2) models with maximal central charge are 
shortened. 

Since A + and iA_, and B + and iB_ are hermitian, the BPS bound is H > \ReZ\. For 
the soliton one has classically that Z is real and H = ReZ, i.e., the bound is saturated. At 
the quantum level, the 1-loop corrections to Re Z are given by expanding ip = tp so \/\/2 + x 
and taking the vacuum expectation values in the soliton vacuum 



Z = Z d + ReJ {u'{ Vbo1 /V2){xx') + ~U"(<p aol /V2){ X xW so i/^} dx 
= Z cl + Rejd x ~U'(tp sol /V2){xx) dx 

Tfl 

= Z d + -Re{{ X {+oo)x{+oo))-{ X {-oo)x{-oo))) (76) 

Since asymptotically (xx) = ( on ly (x*x) is nonzero), there is no correction to the central 
charges. Because there is also no correction to the mass of the soliton, the BPS bound 
remains saturated. 



6 Higher loops 

In this section we repeat the two-loop calculation of the mass of the soliton in the sine- 



Gordon theory |20], [19|], paying close attention this time to possible ambiguities. We begin 
with a review of the method of quantization of collective coordinates, focusing on possible 
ordering ambiguities. For early work on quantization of collective coordinates see |25|, ^6|, 
0, [2]| ^, 0, |3T|, |32|, |33], 0, [3S|. For an introduction see 0. 

To compute the higher loop corrections to the mass of a soliton, one may use standard 
quantum mechanical perturbation theory. One expands the renormalized Hamiltonian into 
a free part and an interaction part, H = H° + Hi nt , and the latter is expanded in terms of 
the dimensionless interaction parameter JhcX/m 2 as H int = J2™=i H^t - 8 This expansion is 



performed both in the soliton sector and in the trivial sector, and the corrections to the mass 
of the soliton can then be evaluated to any given order in hX/m 2 by subtracting the energy 
of the vacuum in the trivial sector from the energy of the vacuum in the soliton sector. For 
the two-loop corrections (themselves of order h 2 X/m since the classical energy M c i of the 
soliton is proportional to m 3 /A) this means that we must evaluate 

= ( S ol|^ sol |sol) - <0|#S itriv iO) 
.( S ol|^«b)(p|^2 iSol |sol) 



p 



E so \ — E p 



8 Terms with n quantum fields rj contain in addition a dimensionless factor (luL/c) n / 2 . We set c — 1 but 
keep h when useful. 
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^/ (OlgStrivlpXplggtnvl ) (77) 

p E — Ep 

Here |sol) is the ground state in the soliton sector (the soliton vacuum) with classical energy 
E Bol = M cl , |0) is the ground state in the trivial (nontopological) sector with vanishing 
classical energy, \p) are the complete sets of eigenstates of H^J and H^ v with positive 
energies E p , and the sums extend over all excitations but do not include the ground state. 
Hence E so \ — E p and Eq — E p never vanish. The fields for the quantum fluctuations in the 
trivial and the topological sectors are expanded into modes with creation and annihilation 
operators, and both |0) and |sol) are annihilated by the annihilation operators. 

6.1 The quantum Hamiltonian 

To apply this approach to the soliton in sine-Gordon theory, we begin by defining the sine- 
Gordon action| 

The action in the trivial sector is obtained by expanding about the trivial vacuum. Since 
for sine-Gordon theory the latter is given by $=0, we obtain = $ + 77 = 77, and 

C=\r ] *-\{r ] y-\mW + ^+... (79) 

In 1+1 dimensional linear sigma models only mass renormalization is needed, = m 2 +Sm 2 , 
and at the one-loop level 5m 2 is fixed by requiring that the graph with a seagull loop and 
two external rj fields cancels the contribution from — (l/2)5m 2 7] 2 . This yields 

X h hX r A dk . . 

Sm 2 = - Y = — / ; 80 

2^2cu vac L An Jo ^k 2 + m 2 V ' 

The mass m is thus the physical mass of the meson at the pole of the propagator to this 
order. 



A complete counterterm which removes all equal-time contractions is [[15 

AH = ( e 5m2 / m2 - 1)^ / + °° (l - cos —Mx)) dx (81) 

A J-00 \ m J 

The corrections to the physical mass at higher loop orders are then finite, and by expanding 
the final result for the soliton mass in terms of the physical mass of the mesons, any ambiguity 
due to defining a AH which differs from ( [Sip by finite terms will be eliminated. In particular, 
the contributions from other renormalization conditions for m, and finite renormalization of 
A and 77, should cancel. 



9 To compare with go) we use their action. It is related to the action in section 5 by the shift 
y/Xtpi/mo — > y/Xipi/mo + 7r. Note that in [po[ mass is renormalized by tuq = m 2 — Sm 2 . 
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The Hamiltonian in the trivial sector is simply 



o-(O) 
^triv 



+00 

dx (82) 



00 



-U 2 (x) + \{r]\x)f + -m 2 r)(x) 2 



^int.triv = -^Ar/ 4 + ... + ^mV + ... (83) 

where Ho{x) is the momentum canonically conjugate to r](x). So H^ triv contains only terms 
for even n and to obtain the contributions of the trivial sector to the two-loop corrections to 
the mass of the soliton we must evaluate — A (0 177 4 1 0) + j5m 2 (0\r] 2 \0). Note that there are 
no ordering ambiguities in the Hamiltonian of the trivial sector. 

To obtain the Hamiltonian in the soliton sector, one must use the formalism for quanti- 
zation of collective coordinates ||. Although the final formulas look somewhat complicated, 
the basic idea is very simple: one expands <f>(x, t) again into a sum of a background field 
(the soliton) and a complete set of small fluctuations about the background field, but in- 
stead of simply writing <f>(x,t) = 4> so \(x) + *E,q m {t)r] m {x) where r) m (x) stands for all modes 
(eigenfunctions of the linearized field equations), one deletes the zero mode for translations 
from the sum, and reintroduces it by replacing x by x — X(t) on the right hand side of the 
expansion of <fi. For small X(t), the expansion of <f) so \{x — X(t)) into a Taylor series gives 
so i — X(t)cj)' sol (x, t) + and since (f>' so i(x, t) is the translational zero mode (the solution of the 
linearized field equations with vanishing energy), one has not lost any degrees of freedom. 
Hence one substitutes 

<f>(x, t) = sol (x - X(t)) + 52'q m (t)Vm(x ~ X(t)) (84) 

into the action in (|78|) , and using the chain rule, one finds an action of the form of a quantum 
mechanical nonlinear sigma model (but with infinitely many degrees of freedom) 

L = u'guiu)^ - V(u); u 1 = {X(t), q m (t)} (85) 

The metric gu is given by 

r d(f>(x, t) d(f>(x, t) . , 

9IJ = J M du' dx m 
and contains space integrals over expressions which depend on q m (t),rj m (x) and 4> so \(x), but 

not on X(t) due to the translational invariance of the integral over x. The Hamiltonian is 
then simply given by 

H = n ig IJ (u)iTj + V(u); Wj = {P(t), jc m (t)} (87) 

where g IJ (u) is the matrix inverse of the metric gu(u) and P(t) is the center of mass mo- 
mentum (the momentum conjugate to while 7T TO (t) are momenta canonically conjugate 
to q m {t). 

Classically, this is the whole result. One may check that the equal-time Poisson brackets 
{Q,P} = l,{g m ,vr n } = 5™ imply {(j)(x),U (y)} = 8{x - y) where U (x,t) = j>(x,t), and 
vice-versa. Hence, the transition from <p(x, t) and n (x, t) to {X(t), q m (t)} and {P(t), n m (t)} 
is a canonical transformation. It is useful to recast the "quantum mechanical" Hamiltonian 
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(89) 



Hi 



U 



one can 



in ( |5?D into a form which resembles more the Hamiltonian of a 1 + 1 dimensional field theory. 
To this purpose we introduce fields constructed from q m and 7r m as follows 

T](x,t) EE Y,'V m (t)Vm(x-X(t)) 

ir(x,t) ee ^2'ir m (t)r) m (x - X(t)) 

By combining the ir m and g m with the functions i] m (x) which appear in g 
write the complete Hamiltonian only in terms of the fields t]{x,t) and tt(x, t) and the back- 
ground field 4> so \(x). To simplify the notation, we introduce an inner product (/, h) = 
f-££ f*(x)h(x)dx. Note that the functions r/ m which parameterize the small fluctuations are 
orthogonal to the zero mode <j)' sol since they correspond to different eigenvalues of the kinetic 
operator 

(0 , sol ,r ?m ) = O (90) 
a result we shall use repeatedly. The classical Hamiltonian density Ti. = T 00 is given by 



U 



^uKx^ + ^'ix^r+v^) 

\-k 2 (x - X(t),t) - n(x - X(t),t) 



[p + M)? 



+ 



2M^l + ( V ',4>' aol )/M d y 



M d [l + W,<lU/M d . 



- x(t)) 



+ -[ v '(x - x(t),t) + tf^ix - x(t))] 2 + v{4> aol + n) 



(91) 



A great simplification occurs in the Hamiltonian H = f_™ Tidx because due to the orthogo- 
nality of the zero mode <j)' sol to the fluctuations r] m , the complicated second term in 7i cancels. 
There should be no terms linear in the fluctuations rj and it and the collective coordinates 
X, P in H (i.e., after integrating 7i over x). That this is indeed the case follows from the 
field equation 

Ci = H0soi) (92) 
The classical energy of the soliton at rest is given by 

"+00 



M, 



cl 



(93) 



which follows from equipartition of energy 

1 

2 1 

Thus we arrive at the following expression for the classical Hamiltonian in the topological 

sector 



ti) 2 



soW 



(94) 



#1°; = M d + f 

J — i 



\v{x,ty + \r ] >{x,ty+ 1 -r?V"{ 



y sol. 



dx 



cl 

int,sol 



A12 



[P + rf)\ 



2M d l + (rf,<fJJ)/M d 



+ 



w + 



> 'sol. 



+ 



dx 



(95) 



(96) 
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All X(t) dependence has disappeared from H due to translational invariance of the integra- 
tion over x. 

We must now discuss the subtle issue of operator ordering in H. We shall consider a 
soliton at rest, so we set P = 0. Furthermore, due to [q m , 7i n ] = %hb™ and (r] m , rj' m ) = 
(since we work in a finite volume, and rj and rf have the same boundary conditions) one has 
the equality (71", rj') = (i]',tt), at least if one considers a finite number of modes in rj and 71. 
However, there are operator ordering ambiguities both in (tt,i]') 2 and also with respect to 
the term (?/, <fi' so] )/M c \ in the denominator. 

In general, one may require that the generators H,P = J T m dx and L = J xT QO dx satisfy 
the Poincare algebra The expressions for these operators are quite complex, and in 

general it seems likely that the operator ordering which leads to closure of the Poincare 
algebra is unique (in quantum gravity, such an ordering has never been found). There is, 
however, an ordering which guarantees closure, and this is the ordering we shall adopt. It 
is obtained by making the canonical transformation at the quantum level. One begins with 
the quantum Hamiltonian in "Cartesian coordinates" (i.e., in terms of the operators Hq(x) 
and 4>(x)). In the Schrodinger representation the operator Tl (x) is represented by d/d<p(x), 
and making the change of coordinates from <p{x) to X and q m , one obtains the Laplacian in 
curved space by applying the chain rule 

1 a r- r ,„ , d 




where a 1 is the set <p(x) and u 1 the set X, q m . If the inner product in a space is given by 



(/, h) = f /*(a)A(a)(n da 1 ), it becomes in u space (/, h) = J f*(a(u))h(a(u))J g^ijldu 1 ). 
With this inner product, the relation between d/du J and the conjugate momenta ttj in the 
Schrodinger representation is not simply ttj = d/du J , but rather 

d 1 = ^Vkrtf^V) (98) 



du 

as one may cheeky Hence the correct quantum Hamiltonian is given by 

1_ 1 r— rr, , 1 

2<? 



H --^7r I ^)g IJ (u)7rj^^+V(u) + AH (99) 



with AH given by (|8l|) . 

It is often useful to rewrite this Hamiltonian such that all expressions are Weyl ordered, 
because then one can use Berezin's theorem and find at once the action to be used in the 



10 The Noether current for the orbital part of the angular momentum J pa is given by = {x p T£ — x a T^) 
so Jox = L = J (x Ti — x\T^)dx. At t=0 this reduces to J xT 00 dx. 

n In Q this derivation of the operator ordering of the Hamiltonian is given, but at the end d/du J is 
replaced by ttj which is incorrect. In E3] the correct quantum Hamiltonian is obtained, but the factors 
g _1 / 4 (u) are produced by "Redefining the Hilbert space so as to eliminate the measure from this scalar 
product...". We claim that the relation (|9^ ) is not a convention or a choice of basis, but is fixed because the 
inner product has been specified. 
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path integral. The result is 

1 



H 
AV 



(ir ig IJ 7ij) w + V{u) + AV + AH 



ft 



djdjg IJ (u) - Ag-^djig^iu^^djg-^iu)} 



(100) 
(101) 



The operator (1/2) (njg IJ r[j^j w is obtained by promoting fl95|) and fl96|) to operators and 

Weyl ordering. Weyl ordering yields then (l/2)((l/4:)n I 7ijg IJ +(l/2)n I g IJ 7ij+(l/4)g IJ ir I 'K j). 
Substituting the expression for g IJ {u) [0 one finds 



i.i 



4g- lf % { 



didjg 



9 1/2 9 IJ djg^} 



dqm dqrt 



(^o, 4L + if 



o 



(^,0') 1/2 dq r ' 







(V o ,0') 7/2 dq" 

where = so i + T], and tpQ = <p' sol / \fM c \ is the normalized zero mode. This leads to 



(102) 



AV 



-2 



(^O,0') 2 

(lPo,V'm)(Vm,Vn){Vn,V') + Oo, V'm) (Vm, V') (Vn, Vn) 



(^o,0O 4 

Further simplifications result by using the identities 



(V>o,0O 



(103) 



(Vn, V) = (Vn, <t>'), (Vm, v'm) = °> J2 Vm(x)Vm(y) = 5(x - y) - 1p (x)lpo{y) (104) 



The final answer for AV reads then 
AV - h -' 



(if 1 2 / Mi£) (^o^o) 

(^O,0') 2 Wo,0') 3 (^,0O 2 



,(iff_ v kVJi 2 



The total Hamiltonian is then the sum of H^} in (|95f) in which no ordering problems are 
present, and (-H^soi) + AV + AH with (-H^soi) given by (|96 ) with the complicated 



(105) 



momentum dependent term Weyl-ordered.This is the result in ||34| . A drastic simplification 
is obtained by rewriting the latter term in a particular non-Weyl-ordered way in such a way 
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that it absorbs all terms in AV except the first one [ 33 1 . This leads to the final form of the 
interaction HamiltonianP] 

^ 2 



H 



1 



int,sol 



8M cl 
h 2 



(P + (i/,7r)), 



1 



+oo 



l + (7/,0' sol )/M cl J 

+ AH 



i)g ol ) 2 dx 



SMZJ-co [! + (//, COW 



+ 



1 



dx 



(106) 



Note that the first term is the square of a Weyl-ordered operator, but is not itself Weyl- 
ordered. 

For the two-loop calculation we are going to perform we set P=0, and we only need the 
terms as far as quartic in 77 and 7t\ This leads to 



H 



int,sol 



2M, 



cl 



+00 



rfudx 



+00 



rj'irdy 



h 2 



8M 2 



+00 



{4>'Ux)fdx 



\/Xm f+°° o, . . \/A . , , , 

+AH — - / rj° {x)sm cp so \(x)dx 

61 J-00 m 



X f + °° 4/ \ VA 

~7T / *7 (^)cos 4> so i(x)dx + 

4! J-00 m 



(107) 



The counterterms are the same in the topological sector as in the trivial sector, but in the 
trivial sector we decomposed <p = $ + 77 = rj, while here we expand <fi = <p so \ + rj. We obtain 
then 



m 



AH = T (e- 



1 



1 — cos —^-d)(x) \ dx 
m 



m 
T 



5m 2 1 ( 5m 2 



f5m \ r Si- cos ^b.^ gfA ^\dx 
\ m z j J-00 I m J 



+ 



5m 2 m r+°° 



T](x) sin ^^4> so \(x)dx 



m 



1 c 2 /" +0 ° 9 / \ v 7 ^ . , . , 
+-om / 77 (a;) cos (p so \(x)dx + 

2 J-00 m 



:io8) 



12 The first term in ( jioq ) can be written for P = as [(tt, r/'/F) + (r//F, tt)] 2 , where F = 1 + {rf , c/)' sol )/M ch 
Weyl-ordering of (tt, rf / F) 2 yields 



*{x)M/F)-L{x) + ^{x){t,'/F,-k)^{x) 



dx + 2{ir,rf/F){rf/F,n) 

Evaluating the difference of the two expressions one needs the following commutators 

V'{y)] = -id y S(x -y) + -^j-<t)' so i{x)<f)" ol (y) 
[tt(x),1/F] = -i4>'Ux)/(M cl F 2 ) 
Straightforward algebra produces then all terms in AV except the first one. 
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The first term is the counterterm for the one-loop graphs and will not contribute to our 
two-loop calculation. 



6.2 The actual two- loop calculation 

Using the explicit expressions for the classical soliton solution 

4m 



so i = -p arctan(e m:c ) (109) 
V A 



one finds M cl = ^- and 



x 



,'VX , \ „ e mx _ e -mx 

sin I so i = -4 T — : , — -ttt; cos so i = 1 



m ' j ( e mx + e -mx^2> m (e mx + e~ mx ) 2 

A 



_ {€o^)fdx = (110) 



8M 2 

Substituting these results into ( 107|J108|) , we find for the Hamiltonian 



„(!),/ v^m /•+- (-4)(e^-e-^0 3 

H(J'i = / ; — 7] {X)dX (HI 

mt,sol g y ^ + e -m*)2 ' W V. V 

{1)jJJ _ £m 2 m (-4)(e--e— ) 

^int.sol " ^ (em , + e - m , )2 VWdX (112) 

<ii = r°°[i- ^ ' , 2 ^ 4 (^ (us) 



flga = 2 5m loo [1 ^ (^r^F ]r? (114) 



OO 



„( 2 ),rest Aft 2 (^m 2 ) 2 4 

^ int ' so1 " 192m + ^A~ m (Ubj 

We now put the system in a box of length L. We expand r](x) into creation and annihi- 
lation operators 



+oo 



V(x) = £ J^(a(q n )E(q n ,x)e ik ^ + h.c.) (117) 

n=— oo V *"^n 

itanhmx + g„ 

^ = + 1 t mh ^]^ (U8) 

q„ = ^;^ = *J + m 2 (119) 
m 

The functions E(q n , x)e lknX ~ tuJnt = t] n (x)e~' luJ " t are eigenfunctions of if^j* and satisfy the 
linearized field equations 

[d 2 x + oj 2 n - m 2 (l - 2 cosh- 2 (x))]r] n = (120) 
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and for large q n they tend to which leads to the familiar normalization factor (2uj n L) 1 ^ 2 
for free fields. At x — ±L/2 we find the phase shift 

Tfl i 

5{k) = 2 arctan{— tanh -mL} (121) 
k 2 

The momenta k n are discretized by adopting antiperiodic boundary conditions^ 

k n L + 5{k n ) = 2im + vr (122) 

As eigenfunctions of a self-adjoint operator on a compact space, the functions E(q n ,x) have 
the same orthogonality properties as plane waves and they are normalized to unity, for 
example 1^1% \E(q n ,x)\ 2 dx = 1. As in the case of the kink 8(k) in ( |121| ) is discontinuous. 



The values of n — . . . , —2, —1, 0,1,2... give solutions of ( |122j ) kL = . . . , — 2ir, 0, 0, 27r, 47r, . . .. 
Again, we see a defect in the n to k mapping: the k = solution is obtained for n = — 1 
and n = 0. If we map k = onto n = — 1, the remaining number n = can be assigned 
to the only discrete solution (the translational zero mode). Then using the Euler-Maclaurin 
formula the sum over n can be converted into integral over k with continuous measure. We 
now evaluate the various terms in (|77|) , adding subsets of terms which combine to cancel 
divergences. 

The i] 4 , 5m? rf, and (5m 2 ) 2 terms in H^ sol and H^ vac 

Since there are 3 ways to contract the 4 7/'s, one finds from the term with yf 



(2) / A f+°° 

( Sol l^intsoll So1 ) = -« 



mx _|_ g— mx\2 



+00/ 

E 



/ J n,p=— 00 

(tanh 2 mx - 1) + (g 2 + 1) (tanh 2 mx - I) + (g 2 + 1) dx 
(1 + aL-^tanh^ (l + g^L-ltanh^ 2o; n 2o; 3 



(123) 



Using the integrals in the appendix, we record the contributions due to terms with none, 
one and two factors tanh x — 1 in 3 separate lines 

(124) 

The first two terms come from the factor 1 — 8/(e mx + e~ mx ) 2 whose x-integral is equal 
to L — — . Clearly, there is a divergence proportional to L which will be cancelled by the 
corresponding contribution from the vacuum sector. The integrals of 1 — 8/(e mx + e~ mx ) 2 
times (tanh 2 x — l) p for p — 1 and p = 2 are given by g 2 - and — respectively, and do 

13 Of course, we differ here from refs. [^9[ who choose periodic boundary conditions. The necessity for 
the antiperiodic boundary conditions has been discussed at length in the preceding sections. 
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not contain divergences which are due to the x integral, but they still contain divergences 
due to the sums over n and p. The terms inside curly brackets are due to expanding the 
denominators (1 + q 2 )L — — tanh ^ ~ (1 + q 2 )L — — . (We have already set tanh ^ = 1 
since the difference vanishes exponentially fast for L — ► oo). Not all these terms vanish for 
L — > oo, but they will cancel with similar terms from other contributions. The sums 



I +oc 1 T 



a; 27r Jo a;(/c) 



1 + !*'(*) 



dk (125) 



will get contributions from S'(k), but by first combining such sums, we shall find that only 
differences like ~ - occur, and this will simplify the analysis significantly. We wrote 
an approximation symbol ~ in (X) -) because the evaluation of such sums has been found 
to be regularization dependent at one-loop level. This is the issue we want to study now at 
the two-loop level. 

The evaluation of (sol|if^^,i|sol) is straightforward and yields 

-8m 2 [ +00 \1 8 1 V ( tanh2mx ~~ X ) + fa* + X ) dx 



2 V— oo ( e m* + e -m*)2J^ ( g 2 + 1)L _A 2uj n 

= ^ 2 (^- 1 )(^E 1 + { J ^E^ + •••}) 

2 m ZL uj mL z 

where the last line contains the contributions from the terms with tanh 2 x — 1. 

The contributions from the trivial sector are only ^Sm 2 {0\r] 2 \0) — ^ (0|?7 4 |0), where rj(x) 
are now simple plane waves and they yield 

(0|iJ triv + AiJ triv |0) = -\5m 2 £ -L _ (127) 
Finally there is the contribution from the term proportional to (5m 2 ) 2 in AH SQ i, it yields 
, ,,l(5m 2 ) 2 f+°° r \/A , . S1 A /V ^lx9 

where we wrote 5m 2 as 

We first demonstrate that the terms due to expanding the denominators (1 + q 2 )L — ^ 
cancel. Consider as an example the last term in the first line of ( |124| ). It cancels the last 
term in the first line of ( |126|) .To see this, write Sm 2 into its original form as a sum over 
modes 

*" 2 -S<^> < 129 > 

One finds for the terms with L~ 2 

= -^(Ei-Ei)(E^) (i3o) 
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Note now that 4 tt is finite, while 



1-1 N 1 AT i 

iKE+E^-E-H 



1 

2^ 



1 



-N 

5(k)_ 
{ L 



)dk -> 



(131) 



is unambiguous and vanishes (the absence of the n = contribution allows the application 
of the Euler-Maclaurin formula as discussed after ( |122| )). Hence, the contributions due to 
expanding the denominator (q 2 + 1)L — — cancel. Similar cancelations occur in other pairs of 
corresponding terms, and we shall therefore only use the terms (q 2 + l)L in the denominators 
of E(x, q n ) from now on. 

The sum of all contributions from the terms with rj 4 , 5m 2 7] 2 , and (5m 2 ) 2 in the topological 
and trivial sectors is then found to combine into differences, except for one term, namely the 
contribution from the first term in the third line in (|124 ) 
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(132) 



Drawn lines represent propagators in the soliton sector and dotted lines propagators in the 
trivial vacuum. In the intermediate expression the first column gives the contributions from 
the 7] A term in H^ sol while the second column gives the contribution from 5m 2 ri 2 in AH so \. 
The third column contains the contributions from the vacuum sector (in the first row) and 
from the term (5m 2 ) 2 in the topological sector (in the second row). We recall that uj and Co 
denote the frequencies in the topological and trivial sectors respectively. We claim that all 
differences cancel. Since we already proved this for (|131|) , we only need to discuss the sum 
t (J2 - — S ^) 2 - Each factor J] - — J2 i is ambiguous but finite, see ( |131| ). Hence the extra 



y factor ensures that also this term vanishes. 



We conclude that all rf and 77 terms contribute only one term 



Am 3 ^ 1 2 
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AQ-K 2 m 



(133) 
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This is the contribution due to the one-vertex two-loop graph in which one only retains in 
both propagators the deviations from the trivial space propagators. The individual other 
terms are ambiguous and divergent but their sum cancels. 

The rf and 5m 2 r] contributions with one intermediate particle 

We next evaluate the contributions to the mass of the soliton which come from the ?f term 
in H^ sol , and the rj term in AH so \. We first take one-particle intermediate states in the 
sums over \p) in (1). Since there are no terms in the trivial sector which are odd in 77, 
these contributions should sum up to a finite result, but we are again interested in possible 
ambiguities. Using 



m\ / 2L\ cosh ±7r<? n 

(see the appendix), and 



. .5m 2 m r+°° -A(e mx - e ~ mx ) , . 

(n\ — =— / — r^-ratx sol) 

y/\ J-00 (e mx + e- mx ) 2 1 ' 

5 m 2 71 ^F n (134) 



1 y> 1 = 2m ( 135 ) 

L V cosh 2 (±7rg n ) vr 2 



we find straightforwardly 



(5m 2 ) 2 

^|(n|r ? term|sol)| 2 /(-cu n ) = ( 136 ) 

Next we evaluate (p|r/ 3 term|sol) with (p\ a one-particle state. This yields a factor of 3 
times an equal-time contraction at the point x. The x-integrals are given in the appendix 
and one finds 

(p, 1 par% 3 term|sol> = V ff (- £ — - -^) (137) 

XP ' 1 " 17 V 2L4m Ycoshi7rg p V L 47rm 2; V ; 

Then it is relatively easy to obtain 

2 

i?e(sol|?7 3 term|p, 1 part)(p, 1 part 1 77 term | sol) 

v ~ u p 

= (5m 2 ) [ l — V - - -?—] (138) 



Finally we evaluate 



|(p, 1 part|r] 3 term|sol)| 2 

= 1 f A y 1 )2, 1 f A y 1 ) A 
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The sum of the contributions from the 5m 2 r] and rf terms with only one-particle inter- 
mediate states is then 

■ ■ + m—T\ + 



— O 

1 -(5m 2 ) 2 + 5m 2 (— Y~- — ) - -L( A y 1 



Am 2mL c<j 67rm Am 4L ' 

1 . A ^ 1 _ A 

'67rmML ^ ur 1207r 2 m 

' ~ 16mL 2 ^ w ~ ^ <^ + 2A-nmL^u ~ ^ 
A 



(140) 

1207r 2 m v ; 

The only nonvanishing contribution from these terms is thus due to the square of the 
finite part of the matrix element of the rf term, and the latter is again obtained by taking the 
deviation from the propagator in the trivial vacuum (the part proportional to tanh 2 mx — 1) 
in the equal-time contraction of two rj fields. All other matrix elements are divergent, and 
none of them contributes. 



The contribution from three intermediate particles and the rjirrjir term 



From ( |111 ) the matrix element with 3 intermediate particles can be written as follows (after 
the substitution e mx = y) 



dy 2 



1 

2m Jo 

(y 2 - 1) 2 



(y 2 



1)4 -i^ 



(y 2 + 1)3 
where Q = q n + q r + q t . Using 

dy 2 y n y lQ -\y 2 + l)-" 
where B is the beta function, one finds 



(y 2 + 1)^ (y 2 + iy 

InQr + QnQt + QrQt) + i 



+ q r + qt) 



y 2 -i 



^q n q r qt 



(y 2 



B( 



(y 2 + 1) 
iQ n + 1 — iQ n + 1 



+ m) 



(141) 



(142) 



1 

2m 
1 

2m 



~oQ + iQ + i;(Q - ^(qnqr + q n qt + q r qt) - Q(q n q r qt) + 



71 



cosh \i\Q 



1 / 2 i 2 i 2\2 I / 2 2 i 22, 2 2\ 

-gK + w r +u t ) + -{u n w r + u) n u) t + u) r u) t ) 



71 



cosh \isQ 



(143) 



Squaring, one finds for the contributions to due to intermediate states with 3 particles 



£ 2 |(?n, ?r, gt|-H"ff,ii|sol) 
n,r,t 

-A 



r(l).J 

int,sc 

dq l dq 2 dq 3 



-1 



2 10 7r6m 7 cosh §7r(gi + <? 2 + 93) 

x [(l+gDH(l+gi) 2 +(l+gg) 2 -2(l+gf)(l+g|)-2(l+gf)(l+gg)-2(l+g|)(l+ 9 l)]2 
(l+g 2 )3/2(l +g 2)3/2(l +g |)3/2 {v ^2 + ^I^| + ^I^f ) 



(144) 
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The prefactor | is needed since we sum over all q n , q r , q t while the 3-particle states are given 
by a) n a\.a\ |sol). 

Next we evaluate the contribution from the / rj'ir J 77V term. Using that 77 = ^ jn {2uj n L)~ 1 
(a n E n exp — iuj n t + h.c.) while n = Y. n {^nL)~ 1 {—ioj n a n E n ex.]i —iu n t + h.c.) one finds 3 
contributions from the 3 possible contractions^ 



(soll^ntiflsol) 



+ / > + 



j ( C n ^ n C r ^ r -\- C r ^_ n C _ n ^ r -\- C n ^— T C — n _ r J (145) 



n a q / j \ ri.ri^ 1.1 1 r. — ii —* — jt.r 1 

64m J ^ V cu n 

where C nir = j(dE n /dx)E r dx. Dashed lines denote n fields. Using that E n and E r are 
orthonormal, it follows that 

11/ \ ITT 

C n>r = imq n 5 nyr + (g n ~ gr) . x. 1 / V ( 146 ) 

Lu n u r v ' smh ^Ti{q n - q r ) 

Since C n)n is odd in n {C n ^ n = imq n + 2iq n / (Luj 2 )) the contribution from C n<n C r ^ r vanishes, 
and in the remainder the term with imq n 5 n)r cancels in the contribution J2 r ,n{~ ~ l)|Cn,r| 2 - 
This leads to 



(sol I J Tj'ir J r/V-termlsol) 



2 10 m 



A / <h,yih n (y/l + ti A fa 2 - glf (147) 



[sinh i7r( 9l - q 2 )] 2 \ yTT^ / i 1 + SiX 1 + d 

According to Verwaest |^0| the sum of these two contributions is 

A 

60iT 2 m 



(148) 



(In ref . [ B3] these contributions were numerically evaluated) . For us the crucial point is that 
both contributions are finite and hence unambiguous. 

Adding ( |133| ), ( 140 ) and ( 148 ) one finds that the sum of these contributions vanishes. 
Hence, the only contribution to the two-loop correction of the soliton mass comes from the 
term -h 2 /(8M 2 ) J((j)") 2 dx in (]10§) and it yields 

M M = (149) 

192m v ; 

All other contributions combine into unambiguous finite integrals which vanish due to factors 
l/L. 

14 Each of these 3 contributions is proportional to ft 2 because each propagator gives a factor Ti, vertices 
do not contribute factors of K, and the energy denominators in (|77j) yield a factor h^ 1 . 
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7 Conclusions 



In this article two new concepts have been introduced in the theory of quantum solitons: 
topological boundary conditions and a physical principle which fixes UV quantum ambigu- 
ities. The main idea underlying these concepts is simple: the problem must be formulated 
before the loop expansion is performed. This means that the mass of the soliton must be 
defined nonperturbatively, rather than as a sum of the classical result, plus one- loop cor- 
rections, plus etc. We define this mass as the difference between the vacuum energies of 
the system with different boundary conditions. It follows immediately that the boundary 
conditions must be formulated for the full quantum field, rather than for small fluctuations. 

The topological boundary conditions are, in fact, better viewed as conditions which put 
the system on a Moebius strip without boundaries. In the literature one usually employs 
periodic boundary conditions in the soliton sector for all the quantum fluctuations (but not 
for the classical field), or a mixture of periodic conditions for bosons and others for fermions. 
All these conditions may distort the system at the boundaries and introduce spurious extra 
energy 0(L°) contributions which obscure the measurement of the mass of a soliton. We even 
found that most boundary conditions in the literature are not compatible with the Majorana 
condition for fermions in the N = (1,1) case. With topological boundary conditions there 
is no spurious energy introduced at the boundaries since there are no boundaries, and one 
obtains the genuine mass of a soliton. 

In the long-standing problem of UV ambiguities of the quantum mass of a soliton, we 
have taken the point of view that this problem should be recognized as the well-known regu- 
larization dependence of loop calculations in quantum field theory. The action is determined 
only up to local counterterms, and one has to introduce renormalization conditions to fix 
those. With solitons, we want to compare vacuum energies in topologically distinct sectors. 
The new principle we have introduced is that if all mass parameters in the theory tend to 
zero, the topological and the trivial sector must have the same vacuum energy (in the infinite 
volume limit). 

This condition follows immediately from a simple dimensional analysis. Indeed, if the 
physical meson mass m and all dimensionful couplings scale according to their dimensions 
as m — > 0, then the mass of the soliton M as a function of these parameters must also 
scale as m. One can take a different, but equivalent, look at the problem: After rotation 
to Euclidean space one finds that the soliton is a string-like interface between two phases in 
two-dimensional classical statistical field theory. Our principle is then related to the well- 
known property of the interface tension: the tension vanishes at the critical (conformal) 



point [pr| . Note again that the formulation of this principle is entirely nonperturbative. 

It is when we do the loop expansion we see that the principle leads to nontrivial con- 
sequences. This is because divergent contributions proportional to the UV cutoff A and 
independent of m may appear (and they do appear in the N = 1 theory with mode number 
cutoff and topological boundary conditions). Our principle can be used as a renormaliza- 
tion condition to eliminate such divergences. Adopting this principle, the calculation of the 
soliton mass becomes unambiguous if one first differentiates the sums over frequencies w.r.t. 
the mass parameter, and then integrates setting the integration constant equal to zero. 



There exist other methods to compute masses in certain 1+1- dimensional models [32 
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35| , |3q , p7| , |38| , |39| , pl| , p2| . For exactly solvable models one can determine the S-matrix 
exactly by assuming that the Yang-Baxter equation holds and that the S-matrix factorizes 
into products of two-particle S-matrices. This program has been extended to the N = 1 susy 
sine-Gordon model |n], |22| and from the result one reads off that to one loop the soliton 
mass is given by 

j ji 

M = M cl (m) - — , (150) 

where m is the physical renormalized meson mass. This is the same result as we have 
found with topological boundary conditions and our renormalization condition at the con- 
formal point. From our point of view, choosing the Yang-Baxter equation (or the ther- 
modynamic Bethe ansatz together with the quantum Lax-pairs approach — the "inverse 
scattering method") amounts to a choice of regularization scheme, which, at least as far as 
the quantum mass of solitons is concerned, appears to be equivalent to our principle. 

However, it may seem that the BPS bound (H) > \ (Z) | is violated at the one-loop level 
because of the negative sign of since the central charge does not receive quantum 

corrections at one loop |§ [l4j (apart from those absorbed by the renormalization of m). 
We point out that it is the unrenormalized expectation value (H) of the Hamiltonian that 
should obey the bound, not the physical soliton mass M, which may (and does in the N = 1 
case) differ from (H) by an m-independent counterterm. Using the mode-number cut-off 
regularization and topological boundary conditions, we have found that 

A in 

<#>mnc = M d (m) + -r-7r ( 151 ) 

4 Z7T 

where A is the ultraviolet cut-off. This means that the bound is observed by (H). This 
happens because N — 1 susy is not enough to eliminate the linear divergent term in (H). 
This term is positive in accordance with the BPS bound. The bound is not saturated, which 
agrees with the observation of Olive and Witten || that the saturation of the bound is 
related to "multiplet shortening" . The latter does not occur in the N = 1 model. 

Therefore, in the N = 1 theory with mode number cutoff we encounter a situation where 
our physical principle M\ m=Q = leads to nontrivial consequences. The subtraction of the 
trivial vacuum energy as in (|J) eliminates bulk volume O(L) contributions (in susy theories 
they are absent anyway). This subtraction cannot, however, eliminate a possible 0(L°) 
m-independent but regularization dependent 0(A) contribution. The required subtraction 
constant, or counterterm, is determined by the condition M| m=0 = 0. Therefore 

M=(E 1 -E )-(E 1 -E ) m=0 (152) 

is a more complete (compared to (^)) definition of the physical soliton mass. It is clear that 
this is exactly the definition implemented by our d/dm calculation 



M 



In = 2 models we do not encounter any of these issues. There it turned out that 
neither the soliton mass nor the central charge receive quantum corrections, hence the BPS 
bound remains intact and saturated. In contrast with the D = 2 N = 1 case where all susy 
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representations, with and without saturation of the bound, are two-dimensional, the BPS 
bound of the N = 2 models is protected by "multiplet shortening" || . 

An alternative UV regularization which one may use is the energy cutoff. In the mode 
cutoff regularization one truncates the sums over the modes. The energy cutoff amounts to 
first converting the sums into integrals over momenta and then truncating these integrals. 



As was shown in [14] this regularization scheme in the sine-Gordon model leads to a result in 



disagreement with the Dashen-Hasslacher-Neveu spectrum [HJ . In the supersymmetric 
sine-Gordon case the energy cutoff would lead to a vanishing one-loop correction (after 
standard renormalization of m), which is in contradiction with existing exact results |^T], 



22J . We examined the two- loop corrections and found that no dependence on the choice 
of regularization appears there. Therefore the difference between the energy cutoff and the 
mode cutoff is purely a one- loop effect. This suggests that, perhaps, a formulation of the 
theory exists where this effect can be described in terms of a quantum topological one-loop 
anomaly. Moreover, the one-loop correction to the mass M does not depend on the coupling 
constant, thus it is, in a certain sense, a geometrical effect. 
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Appendix 



-}dx = — (154) 



) \e mx + e~ mx ) 2 m 

8 2 

1 — -rl(tanh 2 ma; — l)dx = — (155) 

(e mx + e ~ mx ) 2SK ' 3m y ' 

/+oo g A 

\1 - - -rl(tanh 2 mx - ifdx = (156) 

-oo e mx + e - mx ) 2 5m v ' 



/•+oo P mx _ p~mx i , 2 

/ f mXM z^ e im9nX (^mx-iq n )dx = ±*2- (157) 

J-oo (e mx + e mx y 4m costing™ 

r+oo p mx _ -mx 1 

/ t -e im<? " x (tanhmx - iq n )dx = (ql + l) 2 (158) 

J-oo (e mx + e- mx Y V H 1 64m KHn ' cosh \xq n y ' 

(Substitute e mx = y, write iq n e imqnX as ^^ imqnX and put y 2 = | — 1. Then use beta- 
functions and the identities T(l + x) = xT(x) and r(x)T(l — x) — ir/ shi7ra;). 

1 p 1 2m 

L \ (cosh \Kq p ) 2 = ~^ (159) 

Z^(cosh|7rg p )2 = 3^ ^ h e**+T = VT' ( } 

iy < = ^! (lei) 

L p (cosh \nq p ) 2 15tt 2 1 } 
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